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ABSTRACT 

A modelling method having application to petroleum production, soil cleaning, etc. 
for optimizing faster and more realistically the displacement conditions, in a 
porous medium wettable by a first fluid (water for example), of a mixture of fluids 
including this wetting fluid, another, non-wetting fluid (oil for example) and a gas. 
The method comprises experimental determination of the variation curve of the 
capillary pressure in the pores as a function of the saturation in the liquid phases, 
modelling the pores of the porous medium by means of a distribution of capillaries 
with a fractal distribution by considering, in the case of a three-phase water 
(wetting fluid)-oil-gas mixture for example, a stratification of the constituents in the 
pores, with the water in contact with the walls, the gas in the center and the oil 
forming an intercalary layer, determination, from this capillary pressure curve, of 
the fractal dimension values corresponding to a series of given values of the 
saturation in the liquid phase, modelling the hysteresis effects that modify the 
mobile saturations of the fluids effectively displaced in the sample, that vary 
during drainage and imbibition cycles. 



^BACKGROUND OF THE INVENTION 

Field of the Invention 
[0001] The present invention relates to a method for modelling two-phase 

or three-phase flows in a porous medium, in drainage and imbibition. It is based 

on a fractal representation of the porous medium and on an original approach for 

handling phenomena linked with hysteresis (change in the direction of variation of 

the saturations). 

Description of the Prior Art 

1) Experimental studies 

[0002] Experimental determination of the relative permeabilities of a porous 
medium wherein a multiphase fluid flows is not easy. Measuring operations are 
usually simplified by considering that one of the phases is immobile in a state of 
irreducible saturation. 

[0003] The values are for example acquired by means of a well-known 
experimental method referred to as "steady state" for determining relative 
permeabilities, which allows a three-phase fluid to flow with imposed flow rates 
between the phases. The relative permeabilities expressed as a function of the 
two saturations are calculated by applying Darcy's law to each phase. It is not an 
established fact that the relative permeability measurements obtained by means 
of this method are really representative of the fluid displacements and, in any 
case, they take a long time because, at each regime change, one has to wait for a 
state of equilibrium. 
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[0004] Another known method carries out laboratory tests in order to 
determine measurement tables (as shown in Fig.1) relating the relative 
permeabilities and the saturations for each pair of fluids of the three-phase 
mixture. By adjusting experimental production curves, one tries to progressively 
adjust the three-phase relative permeabilities. These data tables are then entered 

into an Athos® type numerical simulator which computes the fluid productions. 

This method is based on the prior acquisition of many experimental 
measurements progressively adjusted by calibration and takes a long time. 

2) Relative permeability models 

[0005] The known empirical model referred to as Stone's model allows, by 
empirical correlations, to predict data relative to a three-phase flow from data 
corresponding to a two-phase flow. It is valid only in case of a high water 
wettability and it is generally considered to be a poor predictor. 

[0006] There are two known types of physical models for modelling three- 
phase flows, based on capillary pressure curves. The capillary pressure curves 
are connected with a saturation (for example that of the mercury injected) and a 
pore radius from which the mercury stops for a given injection pressure, 

determined by Laplace's law, Pinj. = — . 

r 

[0007] A first porous media representation model is described by : 

- Burdine, N.T. :" Relative Permeability Calculations from Pore Size 
Distribution Data", Trans AIME (1953), Vol.198, or by 
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- Corey, A.T. :"The Interrelation between Oil and Gas Relative Pernneabilities", 
Prod. Monthly {^954\ Vol.19, 38. 

[0008] According to this model, the porous medium is represented by a 
bundle of cylindrical capillaries with a radius distribution given by the capillary 
pressure curve obtained by mercury injection. The permeabilities are obtained by 
applying Poiseuille's law to the flow of fluids in these capillaries. 

[0009] This model is based on the representation of the porous medium as 
an assembly of capillaries with different radii. The relation between the volume 
and the radius of the pores is given by the value of the slope of the pseudo- 
plateau. The three fluids are supposed to share the capillaries between them, the 
wetting fluid (water) occupying the smallest ones, the least wetting fluid (gas) the 
largest ones and the third fluid (oil) a zone with intermediate-size pores. It is not 
possible to describe the interactions between the fluids because, in such a model, 
they flow through separate channels. Finally, this model can be useful only if the 
pseudo-plateau covers a wide range of saturations. According to this model, the 
three phases of a three-phase flow move in different capillaries and there is no 
interaction between them. 

[0010] Another known physical porous medium representation model is 

described 

by: 
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- de Gennes, P.G. : "Partial Filling of a Fractal Structure by a Wetting Fluid", 
Physics of Disordered Materials 227-241, New York Plenum Pub. Corp. (1985), 
taken up by 

- Lenormand, R. : "Gravity Assisted Inert Gas Injection : Micromodel 
Experiments and Model based on Fractal Roughness", The European Oil and 
Gas Conference Altavilla Milica, Palermo, Sicily (1990). 

[0011] According to this model, the inner surface of the pores is considered 
isotropic and has a fractal character, and it can be modelled as a "bunch" of 
parallel capillary grooves so that the pores exhibit a fractal cross section. The 
cross section of each pore is constructed according to an iterative process (Fig.1). 
The half-perimeter of a circle of radius Rq is divided into r| parts and each of 
these r| parts is replaced by a semi-circle or groove. At each stage k of the 
process, Nk new semi-circular grooves of radius Rk and of total section Ak are 
created. 

[0012] The fractal dimension DL of the cross section at the end of stage k 
is related to the number of objects Nk generated with the given scale Ik by the 
relation : 

[0013] The fractal dimension can be deduced from a mercury capillary 
pressure c urve a ccording t o t he following p rocedure. M ercury i s i njected i nto a 
porous medium with an injection pressure that increases in stages. Laplace's law 
allows deducing the pore volume, knowing the volume of mercury injected for a 



given capillary pressure and the drainage capillary pressure curve relating the 
injection pressure to the amount of mercury injected and the curve relating the 
proportion of the total volume occupied by the pores and the size of the pores can 
be constructed. In cases where a wetting liquid is drained from the porous 
medium such as water by gas injection, the correlation between the g as-water 
capillary pressure and the saturation of the wetting phase is given by: 
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[0014] The experimental results readily show that the values of the gas- 
water relative permeabilities expressed as a function of the three saturations, 
obtained from the expressions given by the known models and the phase 
distribution modes in the structure of the pores, are far from the measured values 
and therefore that the models concerned prove to be too simplistic to represent 
the complex interactions that take place between the fluid phases. 

[00151 French Patent 2,772,483 (US-6,021 ,662) describes a modelling 

method for optimizing faster and more realistically the flow conditions, in a porous 
medium wettable by a first fluid (water for example), of a mixture of fluids 
including this wetting fluid and at least another fluid (oil and possibly gas). This 
method involves modelling the pores of the porous medium by a distribution of 
capillaries with a fractal distribution considering, in the case of a three-phase 
water (wetting fluid)-oil-gas mixture for example, a stratification of the constituents 
in the pores, with the water in contact with the walls, the gas in the center and the 
oil forming an intercalary layer. It comprises experimental determination of the 
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variation curve of the capillary pressure in the pores as a function of the 
saturation in the liquid phases, from which the fractal dimension values 
corresponding to a series of given values of the saturation in the liquid phase are 
deduced. It also comprises modelling the relative permeabilities directly In a form 
of analytic expressions depending on the various fractal dimension values 
obtained and in accordance with the stratified distribution of the different fluids in 
the pores. A porous medium simulator is used from these relative permeabilities 
to determine the optimum conditions of displacement of the fluids in the porous 
medium. 

[0016] The hysteresis phenomenon relates to the variations in the 
petrophysical properties (relative permeabilities, capillary pressure, resistivity 
index, etc.) observed according to whether measurements are performed during 
drainage or imbibition (these modes respectively correspond to a saturation 
increase and decrease of the non-wetting phase). This phenomenon must 
therefore be taken into account to obtain representative relative permeability 
values. 

[0017] The prior art concerning hysteresis effects in two-phase and three- 
phase media is described for example in the following publications : 

- Land C.S. :"Calculation of Imbibition Relative Permeability for Two and 
Three-Phase Flow from Rock Properties", Trans AIME 1968. Vol.243, 149, 

- Larsen J.A., Skauge A. : "Methodology for Numerical Simulation with 
Cycle-Dependent Relative Permeability", SPEJ, June 1998, and 
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- Carlson F.M. :"Simulation of Relative Permeability Hysteresis to the Non- 
Wetting Phase", SPE 10157, ATCE, San Antonio, Texas, 4-7 Oct. 1981. 

[0018] Fig. 6 typically shows the course of the two-phase permeability 
curves Krw resulting from drainage up to irreducible saturation in wetting fluid (M), 
then imbibition up to residual saturation in non-wetting fluid (NM). The hysteresis 
phenomenon occurs at two levels. At equal saturations Sg, different numerical 
values are obtained and the end point reached is an unknown that depends on 
the cusp point SgM from which the displacement mode is changed. This 
phenomenon is usually attributed to the trapped non-wetting fluid fraction. At 
equal saturations, the same quantity of mobile fluid is not obtained, which distorts 
the flow characteristics. 

[0019] Practically all the models taking account of hysteresis effects involve 
Land's semi-empirical relation : — — = (1) 

where Cl represents Land's constant. This relation relates the initial saturation Sgj 
to the residual saturation Sgr in non-wetting fluid, in order to evaluate the 
saturations in trapped and free non-wetting fluid. Assuming that this relation is 
valid whatever the saturation, it is applied to determine the intermediate mobile 
fractions during displacement In the case of two-phase flows, associating this 
relation with permeability models provides satisfactory results. 

[0020] In the case of three-phase flows, the hysteresis of the relative 
permeabilities Krg takes on a particular form. A displacement hysteresis is 
experimentally observed, as it is the case with two-phase flows (extent of the 
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direction of variation of the saturations), as well as a cycle hysteresis since the 
permeabilities depend on the saturation record. In Fig. 7, the relative permeability 
curves Krg corresponding to a first drainage and imbibition cycle (D1 and 11 
respectively) are distinct from the corresponding curves (D2, 12) of a second 
cycle. 

[0021] The relative permeability model developped by Larsen et al takes 
these two forms of hysteresis into account. Starting from an approach combining 
Stone's model in parallel with Land's formula and Carlson's interpolation method, 
an approach where only the displacement hysteresis is taken into account, 
Larsen et al have introduced an empirical reduction factor that is a function of the 
water saturation, which allows to approximate to the permeability reduction of the 
gas associated with the cycle hysteresis. 

SUMMARY OF THE INVENTION 

[0022] The modelling method according to the invention allows faster 
optimizing and more realistically the displacement conditions, in a porous medium 
wettable by a first fluid, of two-phase or three-phase mixtures including the first 
wetting fluid and at least a second, non-wetting fluid. It therefore provides 
operators with a more reliable tool for evaluating notably the best displacement 
modes of the fluids in the porous medium, in drainage and imbibition. It is based 
on a fractal representation of the porous medium with modelling of the pores by a 
distribution of capillaries with a fractal section, considering a stratified distribution 
of the fluids in the pores, the wetting fluid being distributed in contact with the 



8 



walls and around the second fluid (or around the other two in case of a three- 
phase mixture). 

[0023] The method according to the invention is applicable in many fields 
where fluid flows in porous media are to be modelled in order to optimize the 
conditions of the displacement thereof in drainage and imbibition. Examples of 
fields of application are : 

a) development of an oil reservoir and notably enhanced hydrocarbon 
production by injection of fluids, u sing for example alternate i njections of I iquid 
and gas slugs (a method referred to as WAG). It constitutes an advantageous tool 
allowing reservoir engineers to study also well productivity and injectivity 
problems ; 

b) soil depollution and notably depollution of industrial sites by injection of 
substances such as surfactants in polluted layers; 

c) cleaning of filter cakes by displacement of the substances retained therein; 

d) wood drying ; 

e) optimization of chemical reactions for example by displacement of reaction 
products in a catalyst mass in order to increase the surfaces of contact, etc. 

[0024] The method according to the invention is directly applicable by 
reservoir engineers in order to determine, for example, the most suitable 
enhanced recovery method to be applied to an underground hydrocarbon 
reservoir. The method can also serve within the scope of industrial site 
depollution operations for example. 
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[0025] The method comprises in combination : 

- experimental determination of the variation curve of the capillary pressure 
(Pc) in the pores of a sample of this porous medium in the presence of a wetting 
fluid and of at least one non-wetting fluid (by injection of mercury in a sample 
placed under vacuum for example); 

- determination, from this capillary pressure curve, of the fractal dimension 
values corresponding to a series of given values of the saturation in liquids; 

- modelling the hysteresis effects that modify the mobile saturations of the 
fluids effectively displaced as a function of the number of drainage and imbibition 
cycles undergone by the sample, involving different non-wetting fluid trapping or 
untrapping constants according to whether a drainage stage or an imbibition 
stage is carried out ; 

- modelling the relative permeabilities directly in the form of analytic 
expressions depending on the various fractal dimension values obtained; and 

- entering the relative permeabilities in a porous medium simulator and 
determining, by means of this simulator, the optimum displacement conditions of 
the fluids of the mixture in the porous medium. 

[0026] The method is applied for example for determining displacements of 
fluid mixtures comprising a first wetting fluid, a second, non-wetting fluid and a 
gas, considering a stratified distribution of the fluids in the pores, the wetting fluid 
spreading out in contact with the walls, the gas occupying the center of the pores 
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and the second fluid being distributed in the form of an annular layer in contact 
with both the gas and the first fluid. 

[0027] The method can notably be applied for determining, by means of a 
reservoir simulator, the optimum characteristics of substances added to wetting 
fluid slugs injected in a formation alternately with gas slugs, in order to displace 
hydrocarbons in place, or those of a fluid injected into the soil in order to displace 
polluting substances. 

[0028] Modelling of phenomena by means of the present method has many 
advantages. It allows a better correspondence with the results obtained in the 
laboratory because the physical phenomena are better taken into account. The 
results of the model are therefore better in case of a scale change for example, 
for modelling an application in an operation field. 

[0029] The calculating time is reduced by comparison with the time required 
when tables are used as in the prior methods. Fractal type modelling can better 
deal with the hysteresis effects encountered when using WAG type injection 
processes. 

[0030] The results of the method can furthermore be perfectly integrated 
into many reservoir simulators : 3D, heterogeneous, composition simulators, etc. 

[0031] Exploitation of the results by application softwares is facilitated. It is 
no longer necessary to perform risky interpolations as it is the rule when working 
from discrete values of the result tables in order to draw isoperms for example. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0032] Other features and advantages of the method according to the 
invention will be clear from reading the description hereafter of a non limitative 
example, with reference to the accompanying drawings wherein : 

[0033] Fig. 1 illustrates, in the form of a table, the connections existing for a 
three-phase mixture between the experimentally obtained relative permeability 
values of a fluid and the saturations for two of the three fluids; 

[0034] Fig. 2 is a fractal representation of a pore; 

[0035] Fig. 3 diagrammatically shows the distribution of the phases of a 
three-phase fluid in a fractal pore with the wetting fluid W in contact with the wall, 
the gas phase G spread over the greatest part of the volume of the pore (radius 
Rg)' the oil O being a layer between the wetting fluid and the gas; 

[0036] Fig. 4 shows, as a function of the water saturation, the capillary 
pressure curve of a sample of Vosges sandstone for example, whose local 
slope S is used to determine the fractal dimension of the pores; 

[0037] Fig. 5 shows, as a function of time, the production curves obtained 
experimentally for gas (G), for water (W) and for oil (O), compared with the 
equivalent curves obtained by simulation by means of the method according to 
the invention; 

[0038] Figs. 6a, 6B show the effects of the displacement hysteresis 
affecting the relative permeabilities of the wetting fluid and of the gas, Krw and Krg 
respectively; 
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[0039] Fig. 7 shows the combined effects of the displacement and cycle 
hysteresis observed experimentally on the gas phase; 

[0040] Fig. 8 shows the key curves used for modelling the untrapping 
constant; 

[0041] Fig. 9 shows the result obtained with the model within the scope of a 
WAG type injection; and 

[0042] Fig. 10 shows the validation of the model on a WAG type 
experiment for the recovery curves of the three phases. 

DETAILED DESCRIPTION OF THE INVENTION 

[0043] The method according to the invention allows determination of the 
three-phase relative permeabilities of porous media by using a fractal type model 
of the porous medium, on the basis of an approach described by : 

- Kalaydjian. F.J-M et al. : "Three Phase Flow in Water-wet Porous Media : 
Determination of Gas-oil Relative Permeabilities under Various Spreading 
conditions", 68th Ann. Tech. Conf. and Exh. of the SPE, Houston, Texas, 1993. 

[0044] The method according to the invention comprises, as mentioned 
above, modelling the flow of the phases with distribution of the phases in the 
fractal structure of each pore. In the case of water and oil flowing in a water wet 
porous structure, the oil flows into the volume of the pore. In the case of a three- 
phase flow, there is a stratified distribution, the water, which is the wetting fluid, 
flows along the walls of the pores, the gas circulates in the volume of the pore 
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and the oil flows between the gas and the water. The saturations are calculated 
as the relative surface area in a cross section occupied by each of the fluids. 

[0045] At equilibrium, all the grooves having a radius greater than Rk, which 
is given by Laplace's law Pc = 2 y/Rk. are occupied by the gas, and the smallest 
tubes by the two other fluids (water and oil). The wetting fluid saturation is thus 
expressed as the fraction of the area of the occupied tubes. 

[0046] Calculation of the fraction of the area of the capillaries occupied by 
the water for all the radii between R^ and Roo leads to the following expression : 



and, as Pc = 2 y/R^, the correlation between the capillary pressure and the 
saturation of the wetting phase is given by : 

1 

where Sw is the saturation of the wetting phase. 

[0047] The graphical representation of this correlation in a log-log diagram 
is a straight line starting from the point (Sw. Pc) corresponding to the largest 
capillary of the fractal structure with a radius Rq. One may suppose that : 

- radius Rq, which is first invaded when mercury is injected (Fig.4), corresponds to 
a saturation value of the order of 1/r = 103. Each segment of the capillary 
pressure curve is a part of a line starting from Rq (assumed to be the same for all 
the different segments), corresponding to the aforementioned correlation Pc, Sw 
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Each line has a given slope, a fractal linear dimension can be associated 
therewith. The values of the slope range from 
-1.5 to -3.3 as shown in Fig.4, which leads to values of the fractal linear 
dimension Dl between 1.3 and 1.7; 

- each domain is reached by the mercury for saturations corresponding to the 
place where Rq is found on each line. 

[0048] The saturation of the two liquids when the gas is present in the pore 
is calculated as explained above for a phase : 



assuming that the two liquids occupy the tubes whose radii are less than or equal 
to Rk and the gas, the centre of each pore. The oil saturation is the relative area 
of the cross section of the grooves occupied by oil, whose radius is less than or 
equal to R^. 

Hvsteresis modelling 

[0049] Calculation of the relative permeabilities requires determination, for 
each phase, of the fraction circulating therein and consequently systematic 
estimation of the saturations corresponding to the stagnant fractions. This must 
be done for the two cases studied, for example drainage of water and oil by gas, 
and water imbibition. 

[0050] The original feature of the modelling procedure accounts for the 
hysteresis directly at its origin, that is at the level of the non-wetting phase 
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trapping and untrapping phenomena. In Fig. 8, if the curves could be deduced 
exactly by translation, this would mean that part of the gas has been trapped 
during the secondary drainage and does not take part in the flow. In fact, this is 
not the case but the fact that, at equal gas saturation, the permeability is lower 
during D2 than during 11 means that the trapped gas fraction is larger during 
drainage. In other words, the non-reversibility of the permeability curves can be 
explained by a non-reversibility (hysteresis) between the trapping phenomenon 
and the untrapping phenomenon. 

[0051] Land's formula (Equation 1) is thus kept, but an untrapping constant, 
valid during the drainage stages, different from Land's constant valid during the 
imbibition stages, is introduced. It thus all comes down to modelling the evolution 
of the untrapping constant during the cycles. 

[0052] The characteristics to be taken into account in the formulation are as 
follows. It is considered that Cp = Cl during all the trapping stages because of 
being directly in Land's conditions of application where a single constant is 
enough to describe the phenomenon. When Sgr is low at the end of the drainage 
process, Co is close to Cp. This simply means here that the gas is easier to untrap 
when a small amount thereof is present in the sample, because it can be 
considered that it is present in large pores that can be readily reconnected during 
re-inJection. When Sgr increases, Co becomes greater than Cl, which shows a 
less efficient untrapping process. 



16 



[0053] For high gas saturations, a certain reversibility can be encountered, 
but a trapped gas fraction remains inacessible in small pores. A low single curve 
Krg is thus reached, which corresponds to a pseudo two-phase case where the oil 
is no longer mobile. 

[0054] According to this representation, based on experimental 
observations and working hypotheses, the value of Cd passes through a 
maximum since the untrapping constant is equal to Cl when Sgr is low and when it 
is maximum, when the low-mobility curve is described. The following expression 
allows reconciliation of all the previous aspects : 



Kl-K 



rg min 



K 



D\ 
rg 



(2) 



[0055] This formulation comprises several parameters : 

Cdm : it takes into account the difference between trapping and untrapping, 

Cl :Land's constant, 

Krg min low-mobility curve, 

E : calibration parameter, 

K'rg : value of the relative permeability to gas at the beginning of the previous 
imbibition, 

K'^^rg * value of the relative permeability to gas on the first drainage curve for the 
gas saturation corresponding to K'rg. 



[0056] Whatever the order and the nature of the cycle considered, relation 
(2) and Land's relation (1) allow determination of the trapped and mobile 
saturations by means of Land's formula while taking into account the three-phase 
character of the hysteresis (displacement and cycles) without using empirical 
reduction factors. 



(3) 



and Sg = Sgt + Sgf (4) 

C is equal to Cl or Co according to the displacement mode, 
Sgf : free gas saturation, 
Sgt : trapped gas saturation. 

[0057] Calculation of the relative permeabilities 
A) Relative permeabilities of liquids 

Application of Poiseuille's law to each capillary of the bundle for the phase 
that occupies the capillary allows to calculate the water and oil relative 
permeabilities (Fig. 5). 

[0058] Other experimental studies (Larsen et al) have shown that there is a 
relation allowing to relate the residual oil saturation during the cycles to the 
trapped gas fraction. 

Sor ~ (Sor)sgt=0 ~3Sgt 

(Sor)sgt=o represents the saturation in residual oil left in place in the medium before 
gas is trapped. 
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[0059] If it is considered only the circulating fraction which contributes to the 
hydraulic conductivity, the relative permeabilities for water and oil are expressed 
as follows : 

f^ro = K^oi^PhiiS, -h S^y - {S^ + (1 - a)S^, + {S^r)s,^.,Y] 

R is the reduction factor linked with the trapping of the non-wetting phase. 

[0060] In these expressions, it is useful to mention that : 

- the irreducible water saturation Swi is assumed to be stable, 

- the size range of capillaries occupied by the mobile oil is calculated as the 
difference between the sizes of the capillaries occupied by the two liquids with the 
total liquid saturation Sl = Sq+Sw and those of the capillaries saturated in water 
and stagnant oil, 

- Kro(2ph.) is the value of the oil relative permeability determined by an imbibition 
test with water and oil. When only the water and oil phases are present and since 
the sample tested is water wet, the oil will flow through the section of the pore 
exactly as the gas in a three-phase flow. 

B) Gas relative permeability 

[0061] Since gas is a non-wetting phase, it occupies the central space of 
the pore and it spreads towards the periphery thereof as the gas saturation 
increases, however without coming into contact with the solid wall (Fig.6). It is 
considered that the gas circulates in a single pore whose radius Rg is given by 
the relation : 



Rg = R0+Ri+R2+ +Rk. 
the gas permeability being then given by : 

Kr,=K^^(l-iS,-^S^ry (5) 



with a = . Dl being the linear fractal dimension of the porous medium and 

2-Di 

S|_ the total liquid saturation equal to 1-Sg. 

[0062] The relative permeability model allowing calculating K^v- ^ro 
Krg is installed in a simulator such as ATHOS® or GENESYS©. This allows 
calibration of experiments carried out in the laboratory and optimizing the 
conditions to be satisfied in order to displace petroleum fluids in place in a 
reservoir, either by gas injection or by alternate injection of water and gas slugs (a 
method referred to as WAG), by taking account of the pressure and temperature 
conditions prevailing at the production depth. 

Validation 

[0063] The method according to the invention has been validated by means 
of various types of experiments : 

- gas was injected into porous media containing water and oil, under various 
conditions. It can be seen in Fig. 5 for example that a very good accordance is 
obtained between the production curves of the three phases (water, oil, gas) 
obtained experimentally and those predicted by the reservoir simulator fed with 
the data obtained in accordance with the method, 
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- gas and water were also alternately injected (WAG type injection). Fig. 10 
shows that an excellent accordance is also obtained in this case for the three 
phases throughout the experiment. 
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